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Abstract 

In this note we study the field theory of dynamic isotropic percolation (DIP) with quenched 
randomness that has long range correlations decaying as r"'^. We argue that the quasi static limit 
of this field theory describes the critical point of long range correlated percolation. We perform a 
one loop double RG expansion in e = 6 — d, d the spacial dimension, and (5 = 4 — a and calculate 
both the static exponents and the dynamic exponent corresponding to the long range stable fixed 
point. The results for the static exponents as well as the region of stability for this long range fixed 
point agree with the results from a previous work on the subject that used a different representation 
of the problem [1]. For the special case (5 = e we perform a two loop calculation. We confirm that 
the scaling relation v = ^, is the correlation length critical exponent, is satisfied to two loop 
order. Simulation results for the spreading exponent in d = 3 differ significantly from the value we 
obtain after Pade-Borel resummation was performed on the e expansion result. This is in sharp 
contrast with the result of a two loop e expansion for the spreading exponent for DIP where there 
is a very good agreement with results from simulations for d> 2. 
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I. INTRODUCTION 



For independent site percolation on Z*^, sites are independently assigned to be open with 
probability p or closed with probabihty 1 — p. The subject of percolation is the study of the 
maximally connected sets of open sites, called clusters, where two sites are connected if they 
share a common edge. It could be proved that there exists a critical p, denoted as pc, such 
that for p > Pc with probability 1 there exists an infinite cluster of open sites and for p < p^ 
the probability of such an event is zero. Independent percolation is the simplest example of 
a system with a second order phase transition. The order parameter is the probability that 
a given site belongs to the infinite cluster, this is zero for p < Pc and non zero for p > p^. In 
the case of independent percolation it could be proved that this order parameter is actually 
continuous. For more general information on percolation consult [2, 3] 

In analogy with the usual second order phase transition wc can define various critical ex- 
ponents [2]. The probability that two distant sites belong to the same cluster for p < Pc 
decays exponentially with a characteristic length (correlation length) ^ (p) . Near pc the cor- 
relation length scales as ^(p) ~ |p — Pc|~'^- The expected cluster size S{p) for p < Pc and p 
near Pc scales as S{p) ~ |p — pd"'''- The probability that a given site belongs to the infinite 
cluster Poo (p) for p > Pc and p near pc scales |p — Pel These exponents are not 

all independent, they satisfy the hyper-scaling relation c? = 2^ + ^ [2]. Another exponent 
which is of interest in this note is the spreading exponent. The shortest path L between 
two points on the infinite cluster a distance R apart scales as i?^" , where Zs is the spreading 
exponent [4] . 

Percolation is equivalent to the Q ^ 1 limit of a Q-state Potts model, that is, the critical 
exponents of the Q-state Potts model in the Q ^ 1 limit coincide with those of percolation 
[5]. This could be seen from the partition function of the Q-state Potts model which could 
be written as a sum over "clusters" which in the Q ^ 1 limit correspond to the clusters 
of the percolation problem. Thus, developing a field theory for the Q-state Potts model 
and performing an RG analysis on this field theory allows us to "calculate" some of the 
exponents of the independent percolation problem. Using such an approach, the critical 
exponents u and 7 and thus also /3, using the hyper scaling relation, have been computed in 
terms of an e = 6 — expansion up to third order [6] . 

The statistical properties of the clusters of independent percolation near a percolation thresh- 
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old can also be studied using the so called General Epidemic Process(GEP). GEP is an 
example of an absorbing state phase transition. For this process a "disease" is spreading 
through a media of susceptible individuals. The susceptible media becomes infected with 
rate dependent on the density of the sick and the density of recovered individuals. After a 
brief time interval the sick recover and are immune after that. The recovered individuals, 
sometimes referred to as debris when GEP is used to describe the spread of fire for example, 
stop the spread of the disease locally. The state with zero density of sick individuals is 
absorbing, i.e. the disease can not spontaneously reappear. The statistical properties of the 
debris clusters that are left behind after the disease has been extinguished are described by 
independent percolation [7, 8]. This description of the independent percolation problem 
allows to probe in addition to the exponents u, 7 and (3 also the spreading exponent Zg which 
is connected to the dynamic exponent z of the DIP field theory by the relation Zs — ^ [7, 8] . 
In this note we study a variation of the independent percolation problem. We study the 
effective field theory of a percolation problem in which deciding if a given site is open or 
closed dependents on its surrounding. Such problems arise naturally in statistical mechanics 
models. For example in the Ising model on Z*^ we can declare a site for which the value of 
the spin is —1 to be open and a site for which the value of the spin is 1 to be closed. If we 
are at an infinite temperature then clearly this defines an independent percolation problem. 
If we are at a finite temperature however there will be correlations for the spin values at 
different sites, and thus we have a correlated percolation problem. 

If the correlations of the occupation variables for the percolation problem are governed by 
a finite correlation length, i.e. they are exponential, then the only effect of this correlations 
is to shift the critical density but they don't change the properties of the phase transition 
i.e. they don't change the critical exponents [2]. 

A way to get different critical exponents is to have correlations that decay at large distances 
as power laws, i.e. the model has an infinite correlation length for the occupation variables. 
In the example of the Ising model given above this will correspond to the system being at 
its critical temperature. 

In [9] Weinrib and Halperin argued that the critical exponents of the percolation transition 
should depend only on the decay of the pair correlation G{r) in such systems. In particular 
for G{r) ~ r~" the transition should be in a universality class which depends only on a 
and d, d- the dimensionality of the problem. Their analysis was based on considering the 
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variance of the particle density in a region of volume They found that if a < d these 
correlations are relevant ii au — 2 < 0, where v is the correlation length critical exponent 
corresponding to the pure percolation problem. Weinrib and Halperin argued that systems 
that satisfy the above criteria belong to a new universality class for which the percolation 
correlation length exponent is i^iong — ~- 

In [1] Weinrib used the mapping of the percolation model to the Q ^ 1 limit of a Q-state 
Potts model to construct an effective field theory of the long range correlated percolation 
problem. The long range correlations of the percolation problem were mapped onto a long 
range correlated disorder in the couphngs of the Potts model. To derive an effective field 
theory Weinrib had to resort to the so called "replica trick" and a cumulant expansion. 
Weinrib performed a rcnormalization group analysis of this field theory and using a double 
expansion in e = 6 — ci and S — 4 — a to one loop he obtained results which agreed with the 
results from simple scaling arguments. 

In this note we take a different approach to the problem. We consider a DIP field theory for 
which the critical control parameter r is disordered with a quenched correlated disorder that 
decays for large distances r as r~". The static properties of the clusters that are left behind 
after the agent has been extinguished arc described by long range correlated percolation. 
Performing the averaging over the disorder for this dynamical model is easy, we do not need 
the replica trick [10]. After that we renormalize the resulting field theory using dimensional 
regularization and minimal subtraction. 

The procedure of dimensional regularization and minimal subtraction results in a power 
series for the critical exponents in terms of e, and S when we have double expansion. For the 
independent percolation problem such an expansion even only to two loops after a Pade- 
Borel resummation gives good agreement for d > 3 with the values of the critical exponents 
obtained from simulation [11]. For the minimal distance exponent the agreement is re- 
markable [12]. One might hope, even if it is not realistic, that such an agreement might 
hold for the correlated percolation problem when the power series are in terms of the two 
parameters e and 5. Higher than one loop double expansion however for such models seem 
to be difficult. In this note we consider the special case e = 6. For such a model we perform 
a two loop calculation. 

The results from the two loop calculation do not agree well with simulation results, this 
is not very surprising. The power series in e might not even be resummable, and even if 



4 



it is, the structure of the problem is much more comphcated than the one of independent 
percolation so higher loop calculation might be needed to get comparable result. 
Another RG approach for calculating critical exponents from field theories is the fixed di- 
mensional renormalization, based on Parisi's "massive" scheme [13]. Such an approach for 
independent percolation for d > 3 up to two loops gives quite good agreement with the 
simulation results [11]. The agreement is better than the one based on the e expansion. 
Such an approach was used in [14] to calculate the critical exponents for 0(m)-symmetric 
Ginzburg-Landau- Wilson model in quenched disorder with power law correlations ior d — 3 
and 2 < a < 3. We have tried to do similar calculation for the long range percolation 
problem for d = 3 and a close to 2. However the calculations did not reveal any fixed point 
other than the pure one even after Pade-Borel resumation. The details of this calculation 
are not presented here. 

The paper is organized as follows. In section II we present the field theoretic description of 
the independent DIP model. We also present a sketch of the renormalization procedure and 
how one extracts the critical exponents from the so called RG equations. In section III we 
discuss our generalizations of DIP where we introduce the long range correlated quenched 
randomness. A first order calculation is performed and the new fixed point corresponding to 
long range percolation is identified. Then we present the results of the two loop calculation 
for the special case e — S. In sec IV we present simulation results for the Voter model. In 
the appendix we provide some details of the calculation. 

II. DYNAMIC ISOTROPIC PERCOLATION 

There are two ways to extract an effective field theory functional which after an RG study 
gives the critical exponents of DIP [12]. One approach is to start from the master equation of 
the microscopic dynamics of a specific model in the DIP universality class. Representing this 
in terms of bosonic creation and annihilation operators and using coherent-states one can 
proceed towards field theory [15]. A second approach is to use phenomenological arguments 
and write down an effective Langevin equation obeying all the requirements and symmetries 
of the theory. We can map this equation into a field theory functional [8]. Following this 
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general principle one arrives at the following effective action for DIP 



I[s,s] 



(frdt I s 



at + A(r-v') + y(2^-s) 



(1) 



where S{r\ t) = X dt' s{r, t') is the density of debris at site r, s(r, t) is the density of infected 
individuals, r is the critical control parameter, A is proportional to the recovery rate, and 
finally s is the response field. 

The naive scaling dimensions of the fields and couphngs for the DIP action are as follows 



(d-2)/2 



((i+2)/2 



(6-d)/2 



where is an arbitrary scale of length an time. We see that the upper critical dimension dc 
of the theory is 6, that is for o? > 6 the theory is asymptotically free. 

In the calculation of the Green's functions of a general field theory ultraviolet divergences 
arise, also infrared divergences if we are at the critical point. For a renormalizable field 
theory the divergences can be removed by absorbing them in the bare coupling constants 
and fields. For the DIP field theory we define renormalized fields and couplings as follows 
[12] 

u = G^g'^Z^Z^-^ix-\ \ = {ZZfl'^Z-^\ 



where, e — 6 — d, — The renormalization constants Z — Z..(ti, /x/A, e) can be 

chosen in a UV-renormalizable theory in such a way that 

G'nn({^^ 0, rr, u, Xr, li) = hm Z-^/^i-^/^G({r, t}, r, g, A, A) (2) 

A— >oo 

and G^j^ in (2) are finite and well defined, here G{...) are the Green's functions of the 
theory. We have regularized the field theory in (2) by introducing a high momentum 
cutoff A. For d < dc the critical theory has IR singularities and for d > dc the theory has 
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UV singularities. Indeed the problematic UV and IR singularities are linked precisely at 
d — dc- What is important is that the determination of the Z factors coming from the 
UV divergence provides information of the critical IR singularities and thus on the critical 

exponents [12, 16]. 

In the explicit calculation that we perform we fix e > and take the continuum limit 
A — > oo and we require that the Z factors absorb the e poles. This procedure is called 
minimal subtraction. Note that in such a calculation Tc is set to zero. By requiring that as 
e — > the theory gives finite results we can calculate the exponents as power series in e. 
The bare Green's functions are independent of the renormalization scale /i therefore from 
(2) follows the Renormalization Group (RG) equation: 



■ d_ 
^ djji 



d d (9 1 - - ' 

^Xr— + kTr— + -{NT + NT) 



where 



dXr 



I3{u) = n 



dTr 



du 



N,N 



du 



T{u) = n 



dlnZ 



k{u) = n 



dfj, 

dlnTr 



f{u) = n 



{r,t,Tr,U,Xr,IJ.) ^0 (3) 



dlnZ 



dfjL 



d\n Xr 



dfjL 



(4) 



The partial differential equation (3) can be solved employing the method of characteris- 
tics. After solving the RG equation and employing dimensional scaling one arrives to an 
asymptotic form, long distance , long time of the Green's function from which the critical 
exponents could be derived [12]. The critical exponents of the percolation problem are 
given by 

^ = 2 - f(u*) 



V 



2 - kiu* 



z = 2 + eK) 

where u* is a stable fixed point of the RG, that is j3{u*) = and [u*] > 0. 
An RG study through an e = c? — 6 expansion for DIP results in exponents which agree 
with the exponents obtained from an e expansion of the Q-state Potts model in the Q ^ 1 
limit. The reason why this is so is given in [7]. 
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III. CORRELATED DYNAMIC ISOTROPIC PERCOLATION 



Let us introduce a variation of DIP in which the critical control parameter r that governs 
the strength of the infection is itself position dependent variable r + c(r), with c(r) some 
random field. 

If we take static Gaussian distributed disorder with correlations < c(ri), c(r2) >~ fS{ri — 
and zero average, f the strength of the disorder, and we perform the average, we observe 
that the scaling dimension of / is 4 — d. This is an irrelevant perturbation near 6 dimensions 
so we expect this kind of disorder not to change the critical behavior. 

If however we assume Gaussian disorder with correlations < c(ri)c(r2) >~ / , ^ ,„ and zero 
average, then the scaling dimension of / is 4 — a which is a relevant perturbation for a < 4. 
The explicit form of the functional is 



5* + A (r - v') + Y (2-5 - s) 



[s, s]^ I d'^rdt 

^ J dtij dt2 J d'^n J d'^r2s{ruti)s{ri,ti 



(5) 



-s{r2,t2)s(r2,t2) . 



If we are only interested in time independent quantities (emerging ast ^ oo) it is convenient 
to go to the quasi-static limit [12]. Taking the quasi-static limit amounts to switching 
the fundamental field variable from the agent density s to to the final density of debris 
(f){r) := S{r, oo) = X s{r,t)dt that is ultimately left behind by the epidemic and the 

associated response field 0(r) = s(r, 0) [12]. 

The structure of the action allows us directly to let 



s{r,t) — > 0(r), 0(r) = A 
This results in the quasi-static Hamiltonian: 



s(r, t)dt 



H 



4>,<P 



d^x 



r-V^ + f(0- 



--/ / d^X, / d''x2(l>{xi)cj>{x)- 



\Xi - X2\ 



-(P{X2)4>{X2) 



(6) 



In addition to the rules coming from the Hamiltonian above we have to specify that closed 
propagator loops are not allowed. 

It is more convenient to carry the calculations in momentum space. The Fourier transform 
of the interaction vertex g{x) ~ is g{k) — v + wk°'~^ for small k [9]. As discussed at 



8 



the beginning of this section v is irrelevant and will be ignored, w > 0. We now absorb w 
in the definition of /. 

Using again the arbitrary inverse length scale n by inspection we obtain that 

r 4-a 



The upper critical dimension is 6 and / is relevant for a < 4. The propagator is 

= ^ (7) 

and the vertices are given hy Ui — —U2 — g and V — -^i^. To extract the divergences we 
have only to calculate the one particle irreducible diagrams denoted here by F. Inspection 
of the naive divergence of the one particle irreducible diagrams show that they arise only 
in the diagrams contributing to F^'^, F^'^ = — F^'^ and F^'^. Here the first index is the 

number of amputated external legs and the second is the number of amputated legs. 
The vertex functions are considered as functions of external momenta and we require that 
Fi'2(0), F2'2(0), Fi'i(O) and 



dp 

scheme 



are finite. The model is renormalizable by the following 

p2=0 



u - G,i,-^Z'Z-'g\ 
V = Fsfi-'Z^Z;'f. 

where Fs = d/ir/A 5 = 4 — a. Evaluating to one loop order the divergent diagrams 
and using minimal subtraction and double expansion, where now we require that the Z 
factors absorb both e and 6 poles, we obtain the following results 

Z 1 ^ 

6e 65' 

^ u 2v 

e d 
4u 12v 

Zu = l + ^, 

2u Av 
Z,^l + — . 

e 



This gives us 
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(3{u) — u{—e + -u — lOv), 

5 8 
^^v) = v{-S + 2« - ^v), 

. . 5u — 8v 
k{u,v) = , 

T(«, .) = -- + -. 



A nontrivial fixed point, which corresponds to long range correlated percolation is obtained: 

. 155 - 4e 



u 



11 ' 
^ 215 -lOe 



which finally gives us 
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7^2 <^ 
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One could carry analogously to Weinrib the stability analysis for the different fixed points 
and he will arrive at the same conclusion as in there. Our model however allows us to 
compute the spreading exponent as well. In order to do this calculation we have to go back 
to the dynamical model. Fortunately to extract the dynamical critical exponent we have to 
only calculate F^'^, as a function of the external momentum q and frequency cu From the 



renormalization of the derivative 



1.1, 



we obtain 

=10=0 



4:6 S 

and from this we conclude that 



(ZZ)'/' = 1 + ^ - ^ (9) 



,v ) = -- 777 + ^ = ^77-77' (10) 



12 3 11 44' 

and thus 



z^2 + ^{u*y)^2-^-^^. (11) 
10 



We are interested in obtaining estimates for the critical exponents for long range correlated 
percolation for d> 3. It is quite remarkable that such estimates for independent percolation 
in c? > 3 coming from an e expansion up to two loops agree well with simulation results [11]. 
The agreement for the spreading exponent is quite remarkable [12]. Although it might be 
unrealistic we are curious whether such an agreement might hold for the case of long range 
correlated percolation. Unfortunately it is seems difficult to carry out a two loop double 
expansion in e, 5. We note here that we have performed a fixed dimension renormalization 
for d — 3, and a close to 2, which did not result any fixed point other than the pure one 
even after Pade-Borel resummation was performed. 

We have performed a two loop expansion for the case e = 6, in this case it is just an 
expansion in e. Such models arise naturally when the correlation are expressed in terms of 
the probability that a random walk starting at a given site will hit the origin, this probability 
for d > 3 is proportional to j^p^ , thus a — d — 2. Examples for such models are the Voter 
Model and the Massles Harmonic crystal in d > 3 [17]. Prom the e expansion we obtain 



1/. , ,69 39, 2/ 59 15, 2/ 145 22, 

-{4u - 6v) + uv{- - -) + u\-- + -) + v\-— + 



2, , , 91 32 , 2/ 47 11 , 2/ 91 10 , 

^ \ ,u V . ,71 3 , o, 37 11 , o, 139 5 , 

Z = 1 + -(_-_) + „„(_ - _) + n\-— + ^) + + —^), 

^ . 1/ N / 91 13 ^ 2/ 47 9 , 2/ 91 7 , 



This gives us 



/ 7 ^ 671 2 757 731 2^ 

Pu — (— e H — u — ov u H vu V )u 

^ 2 72 24 36 ' 

5 4 193 2 293 67 2^ 



, 1 1 37 2 71 139 2 

l{u,v) = — M H — v^ u vu ^ V 

^ ' ' 6 3 216 72 108 

, 5 2 193 2 293 67 2 

k(u,v) — -u V u H uv V 

^ ' ^ 6 3 108 72 54 

We identify a long range stable fixed point: 

* I 59 2 

^ =^+88^' 

* 1 131 2 
V — -e-\ e 

2 176 
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,this results in 

k{u*,v*)^^e, (12) 
Prom the dynamical part of the calculation we obtain 

,^-,1 , 1,3 227 2 569 91 n 1 , 5 , , 9 , 

{ZZ)-2 = 1 + -i-u - —u' + —vu -v--v' + -uvlog{2) + - log(2)ii2 - - log{3y) 

1 /51 2 83 7 2^ 
+ ^^(32" -24""+6'' >- ('^' 

This gives us 

7 2 ,1747 5 , 9 , o /427 1, 67 o 

^ __,+ _, + (___ _ iog(2) + - log(3)).^ - (l44 + 4 log(2))^.. + 

For the dynamic exponent, and consequently the spreading exponent, for the long range 
fixed point we finally obtain 



To summarize, the two loop expansion gives for the correlation length critical exponent 
1/ = |. For the ratio of critical exponents ^ we obtain to one loop 2, compare to the result 
of 1.8 in [17], but there is a big correction of — coming from the two loops. From a 
Pade-Borel resummation of the series for Zg we obtain fti 1.6 for d = 3. In the next section 
we report on simulation results for Zg for the Voter model percolation problem on [17]. 



IV. SPREADING EXPONENT FOR THE VOTER MODEL PERCOLATION 

To obtain the spreading exponent for independent percolation one resorts to the so called 
Lcath algorithm. This corresponds to growing a cluster from a single seed [4]. One could 
stop the growth after a certain number of "steps" or after the cluster hits a certain boundary, 
the first is more natural. For the Voter model percolation problem this approach is not pos- 
sible, we can not grow single clusters since the occupation probabilities are not independent. 
We use the algorithm introduced in [17] to simulate the Voter model. We pick a site and 
we decide that it is going to be occupied, this is our seed. Then we run our algorithm for a 
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1 1.5 2 2.5 

log(R) 

FIG. 1: Plot of log(i?) versus log(L) for the 3d voter model at p = 0.1. The slope of the straight 
line gives Zg ~ 1.32 

cube that is centered at that site but in addition to the rules detailed in [17] when a random 
walker hits the center we freeze it and assign all of it ancestors occupied in the percolation 
problem 

We simulate the Voter model at its critical density p — 0.1, the results are presented in 
Fig. 1. In fact not much fluctuation in the results for Zg is observed for p e [0.9, 0.11]. We 
conclude that Zs ~ 1.32. This is smaller than the exponent of independent percolation which 
is Zs ~ 1.37 [4]. This result clearly does not agree with our e expansion result. 

V. CONCLUSION 

We have investigated the field theory of quenched correlated disordered DIP field theory 
with disorder correlations that decays at large distances r as r~". We have identified a long 
range stable fixed point in a one loop double expansion in e = 6 — d and 5 = 4 — a. Our 

results agree with the results obtained in [1] using a different representation of the problem 
as well as different RG scheme. 

For the special case e — 5 we have performed an expansion to 2nd order in e. For the 
correlation length critical exponent we have obtained u — |, for the ratio of critical exponents 
^ the results is 2 — and for the spreading exponent Zs the result is 2;^ = 1 + |e + + 
^ log(2) — ^ log(3))e^. For d — 3 after a Pade-Borel resummation of the series we obtained 
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Zs ~ 1.6. Prom a simulation of the Voter model in 3d we have obtained Zg ~ 1.32. 
We have also performed a fixed dimensional renormalization at d = 3 and a close to 2, the 
details were not reported in this note, and we observed no fixed points other than the pure 
one. It would be interesting to perform such a calculation in the case of c? = 5 and a = 3 to 
see if a stable long range fixed point will appear and how would the results compare with 
the result of the e expansion. 
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VI. APPENDIX 

To one loop in the quasi-static hmit the diagrams that contribute to the different V s are 
hsted below. 



Propagator 



1 



(15) 



Interaction vertices 




9 



(16) 
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We evaluate those using dimensional regularization. 
For the dynamical theory the propagator and vertices are 

Propagator < e{t) exp{-X{T + k'^)t) (17) 



Interaction vertices -^-( — X^g6{t — t) 




AV 

Notice the appearance of a time delocalized vertex [8]. 

For the two loop calculation we consider all topologically different diagrams that can be 
obtained with our vertices and propagator, we discard all diagrams which contain closed 
propagator loops. 

The values of the diagrams that appear could be represented as sums of three types of 
integrals, or their derivatives with respect to the parameters, a,b and c. 



/i(a,6, c) 

l2{a,b,c) - . 

J [a + H 

73(0,6, c) 



a + kf){b + ki){c+{ki + k2y) 

r dfikid'^k2 

J (^+^ 



){h + kl){c + kl + {ki + k2Y) 



(a + kl) {b + kl) {c + kl + kl + {ki + k2f) 
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Only integrals of type Ii appear in the calculation of the quasi static limit, while all types 
of integrals appear in the calculation of the dynamic exponent. The integral Ii and its 
derivatives are evaluated at the point a = b = c = 1, integral I2 and its derivatives are 

evaluated at the point a = b = l,c = 2 and integral and its derivatives are evaluated at 
the point a = b = l,c = 3. The dimensional regularized form of the integrals Ii and I3 can 
be found in the literature [18]. 

h{a, b, c) = + I) (a^- + + c^-) (18) 



/3(a, 6, c) = G^(^(24(«'"^ + b'-') - li^"~'b + b'-^a) - -(a^^ + 6^^) 
1 ..143 9 . A.., , . 1 . .4. 1 



" ^ 1^ + ^ log(^))(«'^ + ^'«) + - ^ log(^))c^(a + 6) 
+ log(l))c(a2 + 6^) + log(l) - i)a6c)) 

For our calculation we only need the first, or higher, derivative of /2(a, &, c) with respect to 
c for that we have obtained: 

dHaM _ ^.(1(1. + 1^.) + l(-ilog(2)a^ (19) 



dc ''e2'8 2 ' 8 

+ + log(2)6c - ^a^ log(a) - log(2)62 - ^ log(2)c2 

- h'' log(6) + + - 1 log(2)a6 + ^ log(2)ac + -^ac - he)) 
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